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GEOPOTENTIAL  HARMONICS  OF  ORDER  15  AND  30  FROM  ANALYSIS 
v  OF  50  ORBIT  DETERMINATIONS  OF  1967-102B 

by 

Doreen  M.  C.  Walker 

SUMMARY 

orbital  parameters  of  the  satellite  1967-102B,  Cosmos  184  rocket,  have 
been  determined  at  50  epochs  from  some  3900  observations.  For  21  of  these  deter¬ 
minations,  Hewitt  camera  observations  were  available. 


The 


During  the  time  of  the  orbit  determinations  the  satellite  passed  through 
15th-order  resonance.  The  variations  in  inclination  and  eccentricity  have  been 
analysed  to  obtain  six  lumped  15th-order  geopotential  harmonics,  with  accuracies 
equivalent  to  between  0.2  and  3.3  cm  in  geoid  height  and  four  lumped  30th-order 
geopotential  harmonics,  with  accuracies  better  than  1  cm  in  geoid  height. 


The  lumped  harmonics  obtained  in  this  Report  have  been  compared  with  those 
from  the  GEM  10B  and  10C  models  and  with  those  from  GRIM3-L1  and  the  RAPP  1981 
models.  A  comparison  has  also  been  made  with  the  lumped  harmonics  obtained  from 
a  new  PGS-3337  model  published  by  NASA,  and  the  agreement  _with  this  model  is  very 
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INTRODUCTION 


The  satellite  1967-102B  is  the  final-stage  rocket  which  put  Cosmos  184  into 
orbit.  It  entered  orbit  in  October  1967  and  initially  had  the  following 
elements1:  inclination  81.20°;  nodal  period  97.27  min;  apogee  and  perigee 
heights,  711  and  544  km  respectively;  and  eccentricity  0.012. 

In  1986,  after  some  18  years  in  orbit,  1967-102B  was  approaching  15th-order 
resonance,  ie  when  the  satellite  makes  15  revolutions  while  the  Earth  spins  once, 
relative  to  the  orbital  plane.  In  this  Report  the  orbit  of  1967-102B  has  been 
determined  from  radar  and  optical  observations  during  1986  and  1987,  using  the 
RAE  orbit  determination  program  PROP  in  the  PROP  6  version2.  The  satellite 
during  this  time  was  placed  high  on  the  priority  list  for  observing.  The 
inclination  and  eccentricity  have  been  analysed  over  the  time  of  the  orbit  deter¬ 
mination  to  evaluate  lumped  harmonics  of  order  15  and  30  in  the  geopotential. 

2  THE  OBSERVATIONS  AND  ORBITS 

2 . 1  The  observations 

The  orbit  of  1967-102B  has  been  determined  at  50  epochs,  from  some  3900 
observations,  between  1986  June  14  and  198-7  November  13.  The  number  of  observa¬ 
tions  used  in  each  orbit  determination  is  given  in  Table  1  and  the  source  of  the 
observations  is  also  indicated. 


There  were  three  groups  of  observation?  available.  Those  listed  first  in 
Table  1  were  the  most  accurate,  being  made  by  the  University  of  Aston's  Hewitt 
cameras  at  the  Royal  Greenwich  Observatory  site  at  Herstmonceux  and  at  the  Siding 
Spring  Observatory  site  in  Australia;  these  observations  usually  have  an  accuracy 
of  a  few  seconds  of  arc.  The  second  consisted  of  visual  observations  made  by 
volunteer  observers  who  used  to  report  to  the  University  of  Aston  but  now  report 
to  the  Royal  Greenwich  Observatory,  Herstmonceux;  these  observations  usually  have 
accuracies  between  2  and  5  min  of  arc  in  good  observing  conditions.  The  final 
and  largest  group  is  made  up  of  Navspasur  observations  kindly  supplied  by  the  US 
Naval  Research  Laboratory,  with  accuracies  of  about  2  min  cf  arc. 
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Table  1 

Sources  of  the  observation  used  in  each  orbit 


Orbit 

No. 

Source 

of  observation 

Total 

Orbit 

No. 

Source 

of  observation 

Total 

Hewitt 

Camera 

Visual 

US 

Navy 

Hewitt 

Camera 

Visual 

US 

Navy 

1 

16  (S) 

77 

93 

26 

6  (S) 

73 

79 

2 

1 

90 

91 

27 

86 

86 

3 

88 

88 

28 

3 

93 

96 

4 

76 

76 

29 

5(H) 

26 

51 

82 

5 

5 

60 

65 

30 

6(H) 

16 

41 

63 

6 

5(H) 

6 

58 

69 

31 

10 

85 

95 

7 

21 

82 

103 

32 

2  (S) 

88 

90 

8 

10(H) 

3 

84 

97 

33 

i2.(s; 

65 

77 

9 

11 

9S 

106 

34 

3 

84 

87 

10 

4 

83 

87 

35 

69 

69 

11 

5 

71 

76 

36 

5 

66 

71 

12 

12  (S) 

38 

SO 

37 

15(H) 

2 

60 

77 

13 

65 

65 

38 

5(H) 

5 

59 

69 

14 

64 

64 

39 

71 

71 

15 

67 

67 

40 

6  (S) 

73 

79 

16 

5 

61 

66 

41 

96 

96 

17 

4  (H)  4  (S) 

18 

58 

84 

42 

1 

72 

73 

18 

11 

61 

72 

43 

1 

91 

92 

19 

6  (S) 

3 

70 

79 

44 

78 

78 

20 

12  (S) 

45 

57 

45 

90 

90 

21 

4  (S) 

52 

56 

46 

11 

55 

66 

22 

17  (S) 

40 

57 

47 

6  (S) 

22 

55 

83 

23 

30  (S) 

40 

70 

48 

18  (S) 

22 

53 

93 

24 

2 

52 

54 

49 

18 

80 

98 

25 

19 

46 

65 

50 

6(S) 

72 

78 

S  =  camera  at  Siding  Spring,  Australia.  H  ■=>  camera  at  Herstmonceux,  UK 

2 . 2  The  orbits 

y 

The  orbits  were  determined. at  approximately  10-day  intervals  with  the  aid 
of  the  RAE  orbit  refinement  program  PROP  in  the  PROP  6  version,  and  the  orbital 
elements  at  each  epoch  are  listed  in  Table  2,  with  the  standard  deviations  below 
each  value.  The  .epoch  for  each  orbit  is  at  00  hours  on  the  day  indicated,  and 
the  PROP  program  fits  the  mean  anomaly  M  by  a  polynomial  of  the  form: 

M  =  Mo  V  Mjt  +  M2t2  +  M3t3  +  M^t4  +  M5t5  ,  (1) 


where  t  is  the  time  measured  from  epoch,  and  the  number  of  M  coefficients 
used  depends  on  the  drag.  For  1967-102B,  which  was  in  a  nearly  circular  orbit 
with  a  perigee  height  of  470  km,  only  Mo  -  M2  were  needed. 
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For  all  50  orbits  listed  in  Table  2,  the  standard  deviations  in  inclination 
were  between  0.0001°  and  0.0014°,  the  average  being  0.0007°,  which  is  equivalent 
to  about  90  m  in  distance.  The  values  of  standard  deviation  in  eccentricity 
varied  between  0.000002  and  0.000015,  and  the  average  value  was  0.000007, 
corresponding  to  50  in  in  distance.  The  value  of  £  ,  which  is  the  parameter 
indicating  the  measure  of  fit  of  the  observations  to  the  orbit,  ranged  between 
0.26  and  0.92,  the  average  value  being  0.51.  The  standard  deviation  for  the 
right  ascension  of  the  node  was  never  more  than  0.001°. 
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Values  of  orbital  parameters  at  50  epochs  with  standard  deviations 
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Tabla  2.  (continued) 
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3  THE  EQUATIONS  FOR  15th-ORDER  RESONANCE 

The  theory  for  general  |3:a  resonance  has  been  given  before,  in  Ref  3, 
where  all  the  parameters  used  are  defined.  The  theoretical  equations  at  15th- 
order  resonance  for  the  variation  of  inclination  and  eccentricity  are  as  follows: 


-54 


di 

dt 


n  (1-e  ) 


sin  i  'a 


/  \15 

(!) 


(15  -  cos  i>?ls,15,7G1S/7,o  {S'1  sin  d>  -  ^  cos  0} 
+  15  (7)  ^ie,15,8G16,8,l  {^15°  Sin(<1>  -  w) 


-1,0 

+  CJ5  cos  (d>  -  (0)  j- 


+  (15  -  2  cos  i)  F16,1S,7G16,7,-1 


r-- 1-2  _-l,2 

x  -!sis  sin (d>  +  o»  +  C15  cos  (d>  +  ca)  J- 


_-l,2 


(2e/2)  q  cos 

+  terms  zn  ■{  - - -  .  (Yd>  -  q(o) 

(lql)l 


(2) 


and,  on  using  equation  (21)  of  Ref  3  after  omission  of  e2  terms  that  are 
negligible. 


f-"(!)  f  flS.!S.7G15,7.0  (S”s‘ 


0,1  _°,1 

sin  d>  -  S15  cos 


d.) 


_1  (f)  ?16.15.8G16.8,1  {\s  »*»<*-  W>  +  C°S "  ®>  } 


+e_1  (f)  ?16.1S,7G16,7,-l{C'2  si“«®  +  ®>  +  5is'  C0S  (<t>  +  ®>  } 


+  terms  zn 


(2/2)  lqle'ql'1 

(iql) : 


2 ,  cos 


{q  -  54(k  +  3q)e  }  sin(Yd>  -  qco) 


(3) 
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The  resonance  angle  <t>  is  given  by 

d>  -  co  +  M  +  15  (Q  -  v)  ,  (4) 

where  v  is  the  sidereal  angle,  and  at  exact  resonance  <5=0.  In 

equations  (2)  and  (3)  only  the  three  terms  with  (y,q)  =  (1,0)  (1,1)  and  (1,-1)  are 

given  explicitly. 

_0,1 

The  quantities  such  as  CJ5  in  equations  (2)  and  (3)  are  lumped  geo¬ 
potential  harmonics  of  order  15,  which  can  be  expressed  as  linear  sums  of 
individual  coefficients  C^a  and  S£n  , 


cq' 


k 


Xq,k- 
t  CiB  and 


_q,k 

s 

m 


(5) 


and  these  individual  coefficients  can  be  evaluated  when  enough  lumped  harmonics 
from  satellites  at  different  inclinations  are  available.  The  summation  for  Z 
in  equation  (5)  is  in  steps  of  2,  beginning  at  the  lowest  Z  ,  which  is  either  15 
or  16  for  the  lumped  harmonics  in  equations  (2)  and  (3) ,  see  Ref  3;  and  for  15th- 
order  resonance  k  =  y  -  q  . 

The  orders  of  magnitude  of  the  terms  in  the  equations  for  di/dt  and 

de/dt  can  be  estimated,  as  the  C.  (or  S.  )  are  expected  to  be  of  order 
-  ,  -q,k 

10”5/82  ,  so  the  value  of  C  (or  S  )  can  be  taken  to  be  of  order 

{£(Q4  x  10  /l2)  }  ,  the  Q  coefficients  being  obtained  for  specific  values 

of  (y,q)  using  the  RAE  computer  program  PROF. 

Estimating  the  orders  of  magnitude  of  the  terms  in  equation  (2)  for 
inclination  suggests  that,  if  just  the  (y,q)  =  (1,0)  and  (2,0)  terms  are  used, 
the  likely  error  is  6%  from  the  neglect  of  the  (1,±1)  terms.  In  equation  (3)  the 
(y,< j)  =  (1,1)  and  (1,-1)  terms  are  dominant,  the  (y,q)  =  (1,0)  term  only 
contributing  about  0.5%;  but  the  (y,q)  =  (2,±1)  terms  may  need  to  be  considered. 

4  ANALYSIS  OF  INCLINATION 

Cosmos  184  rocket  passed  through  exact  15th-order  resonance  on  1987  July  30 
and  its  orbit  has  been  determined  either  side  of  resonance,  with  the  rate  of 
change  of  resonance  angle,  <5,  increasing  from  -2.1  to  +1.2  deg/day.  The 
variation  of  <5  and  <1>  ,  the  resonance  angle,  given  by  equation  (4)  are  shown  in 
Fig  1. 
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In  order  to  analyse  the  changes  in  inclination  due  to  resonance,  all  other 
known  perturbations  must  be  removed.  The  SO  values  of  inclination  in  Table  2 
were  therefore  cleared  of  lunisolar  and  zonal  harmonic  perturbations,  by  using 
the  computer  program  PROD4  with  1-day  integration  steps,  and  the  perturbation  due 
to  the  J2,2  tesseral  harmonic,  which  is  recorded  on  each  PROP  run,  was  also 
removed.  When  the  standard  deviation  for  inclination  was  less  than  0.0005°,  it 
was  increased  to  0.0005°  to  allow  for  the  perturbation  due  to  the  neglected 
effects  of  earth  and  ocean  tides. 

These  revised  values  of  inclination,  with  the  standard  deviations  quoted  in 
Table  2  except  where  they  have  been  increased  to  0.0005°,  were  fitted  with 
equation  (2)  in  integrated  form  using  the  RAE  THROE5  computer  program.  This 
program  removes  the  further  perturbations  due  to  atmospheric  rotation  and  luni¬ 
solar  precession  of  the  Earth's  axis.  The  density  scale  height  H  was  taken  as 
55  km,  appropriate  to  a  height  of  530  km,  1.1  H  above  perigee6,  and  the  atmos¬ 
pheric  rotation  rate,  A  ,  was  taken  as  1.2  rev/day7.  The  values  of  M2  were 
altered  to  mean  values,  M2  ,  by  the  technique  described  in  Ref  8. 

The  first  fitting  by  THROE,  with  <Y,q)  =  (1,0)  (2,0)  gave  £  =  2.9  ,  where 
£  is  the  measure  of  fit  parameter.  On  inspection  of  the  PROP  runs  for  some  of 
the  ill-fitting  values  of  inclination,  it  was  found  that  where  PROP  runs  had  more 
than  one  set  of  Hewitt  camera  observations,  relaxation  of  the  standard  deviations 
of  the  set3  away  from  the  epoch  of  the  PROP  run  gave  better-fitting  values  of 
inclinable*  The  need  for  relaxation  is  justifiable  because  of  limitations  in 
the  accu.  .  the  PROP  orbital  model. 

With  these  relaxations  made,  the  value  of  £  obtained  was  1.9.  A  final 
THROE  run,  with  the  standard  deviations  increased  by  a  factor  of  two  for  three  of 
the  values  of  inclination  to  keep  all  weighted  residuals  less  than  2e  ,  gave 
£  *>  1.64  and  the  values  cf  the  lumped  harmonics  were  as  follows: 


9  0,1 

10  Cls  °  -18.04  ±  0.44 


9_0.2 

10  C3Q  “  0.35  ±  Q.S’' 


9-0,1 

10  SJS  -  -1.62  ±  0.33 


9-0,2 

10  S30  -  8.80  ±  1.08  . 


(5) 


A  further  THROE  run,  with  the  (y, q)  =  ( 1 , ±1 )  terms  added,  was  tried,  but  these 
were  not  well  determined  and  did  not  improve  the  fit. 
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The  values  of  inclination,  cleared  of  all  known  perturbations  except  those 
due  to  resonance,  are  plotted  in  Fig  2.  The  theoretical  curve  derived  from  the 
THROE  fitting  with  (y,q)  =  (1,0)  (2,0),  that  gave  the  values  in  equation  (5),  is 
also  shown  as  a  full  line. 

5  ANALYSIS  OF  ECCENTRICITY 

The  50  values  of  eccentricity  in  Table  2  were  fitted  with  equation  (3)  in 
integrated  form  using  THROE.  The  eccentricity  values  were  cleared  of  zonal 
harmonic  and  lunisolar  perturbations,  and  variations  due  to  air  drag,  within 
THROE  before  the  remaining  variation  due  to  resonance  is  fitted  by  the  program. 

Before  fitting  the  values  of  eccentricity,  the  M2  values  were  altered  to 
mean  values,,  M2  ,  as  for  the  inclination,  and  16  values  of  eccentricity  with 
standard  deviations  less  than  0.000005,  had  them  increased  to  that  value  in  case 
there  was  a  significant  effect  from  the  neglect  of  earth  and  ocean  tides. 

In  some  previous  fittings  of  eccentricity  with  THROE,  it  has  been  found 
that  adjustment  of  the  odd  zonal  harmonics  was  needed.  The  same  method  has  been 
used  here  as  in  Ref  9,  ie  to  add  an  increment  AJ3  to  the  J3  value  used  in  the 
PROP  model. 


With  these  adjustments  made,  the  fitting  with  THROE  was  undertaken.  After 
experimenting  with  pairs  of  coefficients,  the  best  result  was  obtained  with 
(Y,q)  “  (1,1)  (1,-1)  and  (2,1),  giving  8  =>  4.6  .  Other  fy,q)  terms  tried  were 
(2,-1) (1,2)  and  (1,-2)  and  these  resulted  in  undetermined  coefficients.  During 
the  various  tittings  it  was  found  that  a  linear  term  improved  the  fit,  so  the 
value  of  scale  height  was  optimized  to  give  a  better  fit.  The  final  run,  with 
(Y,q)  “  (1,1)  (1,-1)  and  (2,1),  AJ3  optimized  at  -0.237,  H  optimized  at  60  km 
and  standard  deviations  increased  by  a  factor  of  two  on  seven  of  the  values  to 
keep  all  weighted  residuals  less  than  2e  ,  resulted  in  the  following  values  of 
lumped  harmonics,  with  £  »  3.310  . 

9  1,0  9  1,0 

10  C15  -  -41.6  ±  1.9  10  s15  “  "50.8  ±  1.8 


9--1.2 
10  C  IS 


-  -68.3  ±  15.4 


9--1.Z 
105  15 


-90.4  ±  16.3  (6) 


9-1.1 

10  C3Q  -  8.3  ±  1.4 


9-1.1 

10  S30 


-  -12.3+0.8 
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The  values  of  eccentricity,  cleared  of  all  known  perturbations  larger  than  the 
standard  deviations  of  the  value,  except  those  due  to  resonance,  are  plotted  in 
Fig  3.  The  fitted  curve  is  that  given  by  the  THROE  fitting  which  produced  the 
values  in  equations  (6)  . 

6  INCLINATION  AND  ECCENTRICITY  FITTED  SIMULTANEOUSLY 

The  values  of  inclination  and  eccentricity  fitted  separately  by  THROE  can 
be  fitted  simultaneously  using  the  RAE  computer  program  SIMRES  developed  by 
Dr  R.H.  Gooding.  This  program  combines  the  results  from  a  number  of  THROE  runs, 
with  the  same  set  of  (Y,  q)  terms,  and  produces  a  single  set  of  coefficients  to 
fit  the  data.  For  this  SIMRES  fitting,  the  results  from  THROE  runs  with 
(Y ,g)  =  (1,  0)  (2, 0)  (1, 1)  (1, -1)  and  (2,1)  were  used.  The  SIMRES  program,  allows  a 
choice  of  weighting,  so  the  contributing  THROE  runs  can  be  given  more  or  less 
weight  according  to  their  accuracy  of  fit,  which  is  indicated  by  the  value  of 
£  . 

The  THROE  fittings  of  inclination  and  eccentricity  with 
(Y,q)  =  (1, 0)  (2, 0)  (1, 1)  (1, -1)  and  (2,1)  gave  £  =  1.458  and  3.321  respectively. 
For  the  SIMRES  fitting,  therefore,  the  weighting  of  eccentricity  was  down-graded 
by  a  factor  equal  to  the  ratio  of  the  values  of  £  on  the  THROE  fittings,  namely 
2.278  (=3.321/1.458) .  The  values  of  the  lumped  harmonics  given  by  the  SIMRES 
fitting  are: 

9  0,1  9  0,1 

10  C15  =  -19.37  ±  0.45  10  SJ5  =  -2.13  ±  0.34 


9—0,2  9—0,2 
10  C30  =  2.  48  ±0.67  10  S30  “  5.98  ±  1.06 


9—1 ,0  g  1 , 0 

10C.,  =  -41.2  ±2.1  10  S.,  ”  -50.2  ±  1.9  (7) 

1  o  15 


9_  -1  , 2  g-1,2 

10  C1S  =  -70.1  ±16.7  10  S  is  =  -85.0  ±  17.6 


9-1,1  9-1,1 

10  C,„  =  9.0  ±1.5  10  S„.  “  -12.3  ±  0.8  . 

30  jO 
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It  is  not  clear  whether  these  values  are  better  than  those  determined  for 
inclination  and  eccentricity  separately,  equations  (5)  and  (6) :  the  standard 
deviations  are  nearly  the  same.  The  fittings,  if  shown  pictorially  in  Figs  2 
and  3,  would  be  indistinguishable  from  the  curves  given  by  the  THROE  fittings. 

7  EQUATIONS  FOR'  INDIVIDUAL  COEFFICIENTS 

The  lumped  coefficients  in  equation  (7)  can  be  expressed  as  linear  sums  of 
individual  coefficients:  see  equation  (5) .  The  Q  coefficients  in  these 
equations  depend  on  the  ratios3  of  the  eccentricity  functions  G$pq  ,  and  can  be 
evaluated  using  the  RAE  computer  program  PROF.  However,  a  correction  factor  has 
to  be  applied  to  the  values  of  Q  ,  as  it  was  assumed  in  the  PROF  program  that 
the  e2  terms  in  the  function  G^pq  could  be  neglected  (see  Ref  3) .  This  leads 
to  large  errors  in  Q  if  e  is  large;  here,  with  e  approximately  0.01,  the 
correction  is  small,  being  less  than  4%.  The  values  of  the  lumped  harmonics  from 
THROE  and  SIMRES  also  need  a  small  correction,  which  has  been  made  in 
equations  (5)  to  (7) . 

The  equations  given  below  for  the  individual  coefficients  have  been 
terminated  when  the  expected  contribution  from  the  coefficients  permanently  falls 
to  less  than  5%  of  the  largest  contribution. 

The  resulting  ten  equations  are  as  follows: 

41  =  515,1S  +  °-441S7,1S  +  °-173519,lS  "  °-002521,lS  “  0-116623,1S 

"  °-180£23.1S  "  °-2U527(15  -  O'210^  ;  (8) 

_0,1 

S15  :  the  equation  is  the  same  as  (8),  with  S  instead  of  C  ;  (9) 


C30,30  +  °-302C32.3O 


-  0 .2980 


40,30 


0-01l534,30  -°-186S6,30“  0-272638.30 
°-280542,30  -  °-235C"44,30  "  O-176^  f 


.  (10) 

_0,2 

S30  :  the  equation  is  the  same  as  (10)  with  S  instead  of  C  ;  (11) 
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4°  -  Se.15  +  °-910C"18/15  +  °-67l520(15  +  °-384522.15  +  °-106524.15 


°-132C26,15  "  °-308C28,lS  -  °‘415C30.15  "  °-455C32.15 


“  °-434C34,15  “  °-368C36,15  ' 


(12) 


-I/O 

S15  :  the  equation  is  the  same  as  (12)  with  S  instead  of  C  ;  (13) 


C2  =  ^6,15  +  +  2-276520,15  +  2-707522.1S  +  2-939624.1S 


+  2-955C26.1S  +  2-7S9C28,15  +  2-380C3O,15  +  1-863C32,15 


t  1.262C34<1S  +  0.636C36(15  +  0.0400^,.  -  0.480C„#15 


“  °-888C42.1S  '  1-163C«,1S  -  l*299C««.l5  “  1-302C«.1S 


-  1-191C50,15  “  1’001C52,15  ; 


(14) 


_-l,2 

S15  :  the  equation  is  the  same  as  (14)  with  S  instead  of  C  ;  (15) 


C30  -  C31<30  +  0.799C33(30  +  0.494C35(30  +  0.194C37>30  -  0.0540,^ 


“  °-233C41,30  "  °’338C43,30  -  °-374C45,30  "  °-354C47,30 


"  °-293C49.30  -  °-208C51.30  ; 


(16) 


JL.l 

S30  :  the  equation  is  the  same  as  (16)  with  S  instead  of  C  ,*  (17) 
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8  APPROXIMATE  ACCURACY  IN  GEOID  HEIGHT 

Equations  (8)  to  (17)  are  useful  in  allowing  an  approximate  assessment  of 
the  accuracy  of  the  lumped  harmonics  evaluated  here,  to  be  expressed  as  an 
accuracy  in  geoid  height.  If  CT  is  the  error  in  the  lumped  harmonic,  the 
corresponding  geoid  height  error,  CTgeoid  >  is  given  approximately  by  Ror/Q*  , 
where  R  is  the  Earth's  radius  and 


Q*  =  {S  (<£** o/*2)2  f  • 


The  values  of  <7geoid  for  each  lumped  harmonic  in  equation  (7)  are  as  follows: 
lumped  harmonic  Ugeold  (cm)  lumped  harmonic  Ugeoid  (cm) 


0,1 

C15 

0.3 

0,2 

Oi 

CO 

o 

0.4 

1,0 

C1S 

1.0 

-1,2 

515 

3.1 

.1,1 

C30 

0.7 

0,1 

S15 

0.2 

.0,2 

^30 

0.6 

J,0 

S15 

0.9 

-1,2 

S"l5 

3.3 

1,1 

s'3o 

0.4 

9  LUMPED  HARMONICS  FROM  1967-102B  COMPARED  WITH  THOSE  OBTAINED 

FROM  COMPREHENSIVE  MODELS 

It  is  interesting  to  compare  the  lumped  harmonics  obtained  from  1967-102B 
with  those  from  comprehensive  gravity-field  models.  The  comparison  is  made  in 
Tables  3  and  4,  the  values  in  equation  (7)  being  cited.  The  five  models  used  for 
comparison  are:  GEM  10B  and  GEM  IOC10;  the  model  produced  in  1981  by  R.H.  Rapp11; 
GRIM3-L112;  and  a  new  Goddard  Model,  PGS-333713,  a  preliminary  version  of  GEM-T3 . 

The  Goddard  Earth  Model  GEM  10B  extends  to  order  and  degree  36,  and  GEM  IOC 
consists  of  the  GEM  10B  solution  up  to  degree  36,  together  with  some  31000  coef¬ 
ficients  of  order  and  degree  up  to  180,  derived  from  analysis  of  altimeter 
measurements  over  the  oceans.  Equations  (10),  (11)  and  (14) -(17)  require 
coefficients  above  36,  some  as  high  as  degree  52,  but  the  lumped  harmonics  in 
Tables  3  and  4  for  GEM  10B  are  truncated  at  degree  36.  The  expected  accuracy  of 
the  GEM  10B  coefficients  has  been  assumed1,1  to  be  3  x  10"9.  The  GEM  10C  accuracy 
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(above  degree  36)  is  certainly  poorer15  and  is  rather  arbitrarily  taken  as 
5  x  10"9.  The  standard  deviations  for  the  GEM  10B  and  GEM  IOC  values  in  Table  4 
have  been  assessed  using  these  accuracies  for  the  individual  coefficients. 

The  comprehensive  geopotential  model  produced  by  R.H.  Rapp  at  the  Ohio 
State  University  in  1981  gives  the  individual  coefficients  to  order  and  degree 
180,  and  is  derived  from  Seasat  altimeter  data,  terrestrial  gravity  measurements 
and  other  data.  An  accuracy  estimate  for  each  coefficient  is  also  given,  so  the 
standard  deviations  for  the  lumped  harmonics  in  Tables  3  and  4  have  been  assessed 
using  these  accuracies.  The  lumped  harmonics  obtained  by  using  the  individual 
coefficients  from  GRIM3-L1  are  also  given  in  Tables  3  and  4.  This  model  only 
extends  to  degree  and  order  36,  so  all  the  equations  have  to  be  truncated  at  this 
value,  as  with  GEM  10B.  The  standard  deviations  are  calculated  from  the 
accuracies  given  for  the  individual  coefficients.  Finally  the  lumped  harmonics 
obtained  from  a  new  Goddard  solution  PGS-3337,  which  is  complete  in  harmonic 
coefficients  to  degree  and  order  50,  together  with  accuracies  for  each 
coefficient,  are  included  in  Tables  3  and  4. 

Table  3 

Values  of  15th-order  lumped  harmonics  from  1967-102B  and 
comprehensive  geoid  models 


g  0,1 

10  cis 

q  0,1 

10  *15 

9  1,0 
10C15 

g  1,0 

10  S15 

9-"1' 2 

10  C15 

9-'1' 2 
10S15 

1967-102B 

-19.4  ±  0.5 

-2.1  ±  0.3 

-41  ±  2 

-50  ±  2 

-70  ±  17 

-85  ±  18 

GEM  10B 

-22  +  3 

-4  +  3 

-60  ±  5 

-41  ±  5 

-94  +  22 

-74  ±  22 

GEM  10C 

-22  ±  3 

-4  ±  3 

-60  ±  5 

-41  +  5 

-94  ±  26 

-63  ±  26 

GRIM3-L1 

-20  ±  5 

-1  ±  5 

-52  ±  7 

-57  ±  7 

-148  +  27 

-89  ±  27 

RAPP 

-18  ±  4 

-6  ±  4 

-50  ±  4 

-45  ±  3 

-141  ±  19 

-91  ±  19 

PGS  3337 

-17  ±  1 

-5  ±  1 

-51  +  2 

-55  ±  2 

-120  ±  12 

-88  ±  12 
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Table  4 

Values  of  30th-order  lumped  harmonics  from  1967-102B  and 
comprehensive  geoid  models 


9  0,2 

10  630 

9  0,2 

10  S_30 

9-1.1 

10C30 

9.1.1 

10S3O 

1967-102B 

2.5  +  0.7 

6.0  ±  1.1 

9  +  2 

-12  +  1 

GEM  10B 

-5  +  3 

11  ±  3 

0  ±  4 

-19  ±  4 

GEM  10C 

-7  +  4 

13  ±  4 

1  +  6 

-14  +  6 

GRIM3-L1 

-1  ±  3 

6  +  3 

-4  ±  4 

-21  ±  4 

RAPP 

-5  ±  2 

7  ±  2 

-4  ±  2 

-13  ±  3 

PGS-3337 

-2  ±  2 

4  +  2 

-3  +  3 

-15  +  3 

.0,1  .0,1 

In  Table  3  the  values  of  C1S  and  S15  from  1967-102B  are  nominally 
much  more  accurate  than  those  from  the  models:  the  agreement  is  good,  with  eight 
of  the  ten  values  differing  from  the  1967-102B  values  by  less  than  the  sum  of 
their  standard  deviations.  The  mean  difference  for  all  30  model  values  in 
Table  3  from  those  of  1967-102B  is  1.1  times  the  sum  of  the  standard  deviations. 

In  Table  4  the  values  for  GEM  10B  and  GRIM3-L1  should  not  be  taken  too 
seriously,  as  both  models  only  go  to  degree  and  order  36,  whereas  equations  (10) 
and  (11)  require  coefficients  to  degree  46  and  equations  (16)  and  (17)  coef¬ 
ficients  to  degree  51.  The  values  from  the  other  3  models  agree  on  average  to 

within  1.4  times  the  sum  of  the  standard  deviations,  with  those  from  1967-102B. 

.0,2  .0,2 

The  values  of  C3Q  and  S3Q  obtained  here  have  been  used  to  determine  an 
improved  set  of  values  of  individual  harmonic  coefficients  of  order  30  and  even 
degree:  see  Appendix. 

10  CONCLUSIONS 

The  orbit  of  1967-102B  has  been  determined  at  50  epochs  from  some  3900 
observations,  between  1986  June  and  1987  November  while  the  satellite  was  passing 
through  15th-order  resonance.  The  average  accuracy  of  the  inclination  and 
eccentricity  for  all  50  epochs  was  equivalent  to  90  m  and  50  m  in  distance 
respectively. 

The  variations  in  inclination  and  eccentricity  have  been  analysed;  six 
15th-order  and  four  30th-order  lumped  harmonics  have  been  evaluated.  These  ten 
values  of  lumped  harmonics  have  standard  deviations  corresponding  to  between 
0.2  cm  and  3.3  cm  in  geoid  height. 
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The  lumped  harmonics  have  been  compared  with  those  obtained  from  five 
comprehensive  models  of  the  gravity  field,  and  have  been  used  in  a  new  evaluation 
of  individual  harmonic  coefficients  of  order  30  and  even  degree:  see  Appendix. 
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Appendix 

NEW  VALDES  FOR  INDIVIDUAL-  HARMONIC  COEFFICIENTS  OF  ORDER 
30  AND  EVEN  DEGREE,  OBTAINED  WITH  THE  AID  OF 
THE  RESULTS  FROM  1967-I02B 

When  values  of  lumped  harmonics  of  a  particular  order  are  determined  from 
resonant  orbits  having  a  wide  variety  of  orbital  inclinations,  it  is  possible  to 
solve  for  the  individual  harmonic  coefficients  of  that  order.  Such  solutions 
obtained  previously  for  order  30  and  even  degree16  have  used  accurate  results 
from  several  satellites  at  inclinations  between  50°  and  74°,  and  also  one  at  90°; 
but  the  only  previous  results  from  an  orbit  at  an  inclination  between  74°  and  90° 
were  at  80°  from  Ariel  3,  a  satellite  of  relatively  high  drag.  The  new  results 
from  1967-102B,  at  81°  inclination,  are  much  more  accurate  (nominally  by  a  factor 
of  about  5)  and  should  lead  to  an  improved  solution  for  individual  harmonics  of 
order  30  and  even  degree. 

The  previous  solutions,  for  six  individual  C (or  S)  coefficients,  of  degree 
30,  32,  34,.. 40,  were  derived  from  11  satellite  equations  of  the  form: 


C30,30  +  °32C32,30  +  Q34C34,30  +  "  °40C40,30 


_0.2 

~30 


±  a 


(A-l) 


where  the  numerical  values  of  the  Q  coefficients  depend  on  the  inclination,  and 

.0,2 

C3Q  is  the  lumped  harmonic  determined  from  analysis  of  the  resonance,  with 
standard  deviation  O  .  These  satellite  equations  were  supplemented  by  six  con¬ 
straint  equations  of  the  form: 


C 


it  30 


(A-2) 


which  express  the  expectation  that  most  of  the  30  values  will  be  less  than 
10 "5/f2  .  There  were  similar  equations  for  the  S  coefficients. 

If  1967-102B  is  to  be  added  to  these  equations,  the  standard  deviations  of 
the  lumped  harmonics  need  to  be  increased  to  allow  for  the  neglected  coefficients 
of  degree  42,  44,  46,  ...  This  was  done  on  the  assumption  that,  for  l  £  42, 

30  "  0.5  x  10'5/22  ,  and  similarly  for  30  ,  giving  a  total  extra  contribu¬ 
tion  to  a  of  1.1  x  10"9.  When  this  is  added  to  the  existing  errors, 
equations  (10)  and  (11),  truncated  at  degree  40,  give 
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<  ti1- 


C  ‘  +  0.302C  -  O.'OllC,.  -  0-.  186c  , 

30,30;  .-32,30  34,30  36,30 


-  0.272C  -  0 .298C  =  (0.35  +  l.*29)  x  10~9  (A-3) 

38,30  40,30 


and  similarly  for  S  ,  with  the  term  on  the  rhs  as:  (8.80  ±  1.54)  x  10”9.  note 
that  the  harmonics  of  equation  (5)  are  used  here. 

The  new  values  for  C(or  S)  coefficients  are  obtained  from  12  satellite 
equations  -  the  previous  11  and  equations  (A-3)  -  and  the  six  constraint 
equations  (A-2) .  As  before,  the  constraints  were  relaxed  when  necessary  to 
ensure  that  the  weighted  residuals  did  not  exceed  1.0,  and  the  standard  devia¬ 
tions  of  the  lumped  harmonics  were  doubled  for  any  satellite  for  which  the 
weighted  residual  exceeded  1.4.  For  the  S  coefficients,  no  relaxations  or 

doublings  were  required,  except  for  those  previously  made.  For  the  C-coef ficient 

_0,2 

solution  the  standard  deviation  of  C3Q  for  Ariel  3  and  for  1967-102B  had  to  be 
doubled,  and  two  of  the  constraint  equations  had  to  be  slightly  relaxed.  Vhe  new 
values  of  the  individual  harmonic  coefficients  obtained,  with  their  standard 
deviations  are  given  in  Table  A.l. 

Table  A.l 

New  values  of  even-degree  3(J  and  3Q 


t 

1<)9£8.30 

10  £8,30 

30 

-1.8  ±  0.8 

7.6  +  0.7 

32 

-7.6  ±  1.6 

4.7  ±  1.6 

34 

-16.4  ±  1.9 

-5.9  ±  1.8 

36 

-11.2  ±  2.6 

4.8  ±  2.3 

38 

-0.9  ±  2.5 

3.0  ±  2.1 

40 

0.0  ±  2.2 

-4.5  ±  2.1 

The  standard  deviations  of  the  values  in  this  new  solution  are  on  average 
25%  lower  than  before,  so  the  addition  of  1967-102B  does  have  a  beneficial 
effect.  The  average  standard  deviation  is  1.85  x  10-9,  equivalent  to  1.2  cm  in 
geoid  height. 
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The  values  of  the  C  coefficients  differ  appreciably  from  the  previous 
set16:  two  of  them  (for  t  -  36  and  38)  have  changed  by  more  than  the  sum  of  the 
standard  deviations.  The  S  values  are  not  significantly  altered. 

The  weighted  residuals  for  the  12  satellite  equations  and  six  constraint 
equations  are  given  in  Table  A. 2,  where  the  symbol  R  denotes  a  relaxation  of 
the  10"s/£2  constraint.  The  overall  measure  of  fit  £  was  0.82  for  C  and 
0.85  for  S  . 


Table  A. 2 
Weighted  residuals 


Satellite  equations 

Constraint  equations 

0,2 

_o,2 

deg 

Satellite 

C30 

S30 

i 

C£,30 

S£,30 

50.6 

74-34A 

0.13 

0.46 

30 

0.16 

-0.68 

56.1 

68-70A 

-0.53 

-1.26 

32 

0.78 

-0.48 

58.2 

63-24B 

0.72 

0.30 

34 

1.00R 

0.68 

65.0 

65-14A 

0.14 

-0.33 

36 

1.00R 

-0.62 

65.8 

71-10B 

-0.49 

0.98 

38 

0.12 

-0.44 

74.0 

70-111A 

-0.38 

-0.31 

40 

-0.01 

0.72 

74.0 

71-13B 

1.09 

0.16 

80.2 

67-42A 

-1.34 

-1.20 

81.2 

67-102B 

0.76 

0.08 

90.2 

71-54A 

-0.13 

0.07 

98.7 

64-52B (H) 

-0.57 

0.81 

98.7 

64-52B(B) 

0.56 

1.14 

Table  A. 3  gives  a  comparison  of  the  new  values  with  the  corresponding 
values  in  four  comprehensive  geopotential  models,  GEM  10B10,  PGS-333713, 
GRIM3-L112  and  Rapp' 3  1981  model11.  The  values  from  these  models  are  believed  to 
be  completely  independent  of  those  from  resonances,  apart  from  Rapp' 3  value  for 
i  -  30,  which  utilizes  some  resonant  results.  The  errors  of  GEM  10B,  PGS-3337, 
GRIM3-L1  and  Rapp(1981)  are  assessed  as  3,2,3  and  2.5  X  10-9  respectively.  As 
before16,  the  comparison  is  terminated  at  l  -  36  ,  which  is  as  far  as  GEM  10B 
and  GRIM3-L1  go. 
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Table  A.  3 

Comparison  of  new  •  30th-order  values  with 
comprehensive  models,  for  -8  <  3  6 


io9 

^£,30 

109£ 

,  30 

■e 

GEM 

10B 

PGS- 

3337 

GRIM3- 

L1 

Rapp 

New  values 

GEM 

10B 

PGS- 

3337 

GRIM3- 

L1 

Rapp 

New  values 

30 

'  -5.2 

■m 

-0.6 

(-3.3) 

-1.8  ±  0.8 

11.1 

5.4 

n 

(7.5) 

7.6  ±  0.7 

32 

-0.6 

-5.7 

-6.9 

-6.7 

-7.6  +  1.6 

Em 

-0.1 

0.5 

4.7  +  1.6 

34 

-11.9 

-18.6 

-23.0 

-22.9 

-16.4  +  1.9 

-3.1 

0.7 

-0.6 

-5.9  +  1.8 

36 

-3.9 

-8.2 

*  -7.6 

-6.0 

-11.2  ±  2.6 

B 

3.4 

6.4 

•  4.8 

4.8  +  2.3 

The  new  values  agree  well  with  PGS-3337 :  the  mean  difference  between  them 
is  2.5  x  10"9  and  seven  of  the  eight  values  are  within  the  sum-  of  their  standard 
deviations.  The  mean  difference  between  the  new  values  in  Table  A. 3  and  the 
corresponding  values  in  GEM  10B,  GRIM-L1  and  Rapp  (1981)  are  5.4,  3.3  and 
3.7  x  10-9  respectively  (excluding  i  =  30  for  Rapp).  As  the  average  standard 
deviation  of  the  new  values  in  Table  A.  3  is  1.7  X  10"9,  the  differences  are 
consistent  with  the  error  estimates  already  quoted  for  the  models.  All  five  sets 
of  values  in  Table  A. 3  agree  in  having  C34  30  as  the  numerically  largest  value. 

The  new  values  should  be  better  than  those  previously  obtained  from 
resonances,  but  they  cannot  be  regarded  as  completely  reliable  because  of  the 
sparsity  of  data  at  inclinations  between  74°  and  90®. 

Until  now  there  have  been  considerable  differences  in  the  values  of  30th- 
order  coefficients  from  different  sources:  the  agreement  between  the  new  values 
from  resonance  and  the  completely  independent  PGS-3337  model  is  encouraging. 

Author's  note 

This  Appendix  has  been  written  in  collaboration  with  Dr  D.G.  King-Hele  and 
a  version  of  it  has  been  published  in  Planetary  and  Space  Science  (Vol  38, 
pp  407-409  (1990)).  The  results  given  here  include  NASA's  latest  PGS  model  for 
comparison  and  this  gives  much  better  agreement  than  GEM-T2,  which  was  used  in 
the  published  version. 
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Fig  3  Values  of  eccentricity  near  15th-order  resonance,  with  fitted  theoretical  curve 
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analysed  to  obtain  six  lumped  15th-order  geopotential  harmonics,  with  accuracies 
equivalent  to  between  0.2  and  3.3  cm  in  geoid  height  and  four  lumped  30th-order 
geopotential  harmonics,  with  accuracies  better  than  1  cm  in  geoid  height. 

The  lumped  harmonics  obtained  in  this  Report  have  been  compared  with  those 
from  the  GEM  10B  and  10C  models  and' with  those  from  GRIM3-L1  and  the  RAPP  1981 
models.  A  comparison  has  also  been  made  with  the  lumped  harmonics  obtained  from 
g  j  a  new  PGS-3337  model  published  by  NASA,  and  the  agreement  with  this  model  is  very 
"  [good. _ 
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